Cosmology in Six Dimensions* 
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We discuss cosmological models in six-dimensional spacetime. For codimension-1 
branes, we consider a (4+1) braneworld model and discuss its cosmological evolution. 
For codimension-2 branes, we consider an infinitely thin conical braneword model 
in the presence of an induced gravity term on the brane and a Gauss-Bonnet term 
in the bulk. We discuss the cosmological evolution of isotropic and anisotropic 
matter on the brane. We also briefly discuss cosmological models in six-dimensional 
supergravity. 



1. INTRODUCTION 

Recently, there have been many observational and the- 
oretical motivations for the study of theories with extra 
spacetime dimensions and in particular the braneworld 
scenario. From the observational side, the current 
paradigm, supported by many recent observations like 
the cosmic microwave background anisotropies [1], large 
scale galaxy surveys [2] and type IA supernovae [3, 4] 
suggest that most of the energy content of our universe 
is in the form of dark matter and dark energy. Although 
there have been many plausible explanations for these 
dark components, it is challenging to try to explain these 
exotic ingredients of the universe using alternative grav- 
ity theories as such of the braneworlds. From the the- 
oretical side, such extra-dimensional braneworld models 
are ubiquitous in theories like string or M-theory. Since 
these theories claim to give us a fundamental description 
of nature, it is important to study what kind of grav- 
ity dynamics they predict. The hope is to propose such 
modified gravity theories, which share many common fea- 
tures with general relativity, but at the same time give 
alternative non-conventional cosmology. 

The essence of the braneworld scenario is that the 
Standard Model, with its matter and gauge interactions, 
is localized on a three-dimensional hypersurface (called 
brane) in a higher-dimensional spacetime. Gravity prop- 
agates in all spacetime (called bulk) and thus connects 
the Standard Model sector with the internal space dy- 
namics. This idea, although quite old [5], gained mo- 
mentum the last years [6, 7] because of its connection 
with string theory (for a review on braneworld dynamics 
see [8]). 



Cosmology in theories with branes embedded in ex- 
tra dimensions has been the subject of intense investi- 
gation during the last years. The most detailed analysis 
has been done for braneworld models in five-dimensional 
space [9]. The effect of the extra dimension can mod- 
ify the cosmological evolution, depending on the model, 
both at early and late times. The cosmology of this and 
other related models with one transverse to the brane 
extra dimension (codimension-1 brane models) is well 
understood (for a review see [10]). In the cosmological 
generalization of [7], the early times (high energy limit) 
cosmological evolution is modified by the square of the 
matter density on the brane, while the bulk leaves its im- 
prints on the brane by the "dark radiation" term [9, 11]. 
The presence of a bulk cosmological constant in [7] gives 
conventional cosmology at late times (low energy limit) 
[11]. The early time modification for example, can be 
interesting phenomcnologically because may require less 
fine-tuned inflationary parameters [12]. 

In the above models, there are strong theoretical argu- 
ments for including in the gravitational action extra cur- 
vature terms apart from the higher dimensional Einstein- 
Hilbert term. The localized matter fields on the brane, 
which couple to bulk gravitons, can generate via quan- 
tum loops a localized four-dimensional kinetic term for 
gravitons [13]. The latter comes in the gravitational ac- 
tion as a four-dimensional scalar curvature term local- 
ized at the position of the brane (induced gravity) [14] 1 . 
In addition, curvature square terms in the bulk, in the 
Gauss-Bonnet combination, give the most general action 
with second-order field equations in five dimensions [16]. 
This correction is also motivated by string theory, where 
the Gauss-Bonnet term corresponds to the leading order 
quantum correction to gravity, and its presence guaran- 
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There has been a lot of discussion about the potential problem 
of the extra polarization states of the massive gravitons regard- 
ing phenomenology (discontinuity problem), but the particular 
model seems to be consistent in a non-trivial way [15]. 



-2- 



tees a ghost-free action [17]. Let us note, however, that 
if the curvature squared terms are to play an important 
role in the low energy dynamics, one may face a difficulty 
when interpreting gravity as an effective field theory (in 
the sense that even higher dimensional operators would 
seem to be also relevant). 

If curvature corrections are included in [7] , the presence 
of the 3-brane gives similar modifications to the standard 
cosmology. In the case of a pure Gauss-Bonnet term in 
the bulk, the early time cosmology is modified by a term 
proportional to the matter density on the brane to the 
power two thirds [18]. If an induced gravity term is in- 
cluded, the conventional cosmology is modified by the 
square root of the matter density on the brane at low 
energies [19]. Finally, if both curvature corrections are 
present, at early time the effect of both the Gauss-Bonnet 
and the induced gravity terms make the universe to start 
with a finite temperature [20] while at late times the ef- 
fect of induced gravity make the universe to accelerate 
without the need of dark energy [21]. 

Braneworld models can also be extended in higher than 
five-dimensions. Wc can consider a (n + 1) brane em- 
bedded in (n + 2) spacetime or a (3 + 1) brane em- 
bedded in (n + 1) spacetime. In this talk we will dis- 
cuss (4 + 1) cosmological braneworld models embedded 
in six-dimensional spacetime (codimcnsion-1 models) and 
(3+1) cosmological braneworld models embedded in six- 
dimensional spacetime (codimension-2 models). 

Six or higher-dimensional braneworld models of 
codimcnsion-1 are considered as generalizations of the 
Randall-Sundrum model. In [22] static and non-static 
solutions in a six-dimensional bulk as well as the stabil- 
ity of the radion field were discussed, while in [23] the 
evolution of the extra dimensions transverse to the brane 
in a Kasner-like metric was studied. 

In braneworld scenarios, contrary to the Kaluza-Klcin 
theories, the extra dimensions can be large if the geome- 
try is non trivial. If the hierarchy problem is addressed in 
the braneworld scenarios for example, the extra dimen- 
sions should be large [6] . In a cosmological context how- 
ever, these extra dimensions are observationally much 
smaller than the size of our perceived universe. Initially 
the universe could have started with the sizes of all di- 
mensions at the Planck length. Then, a successful cos- 
mological model should accommodate in a natural way 
a mechanism by which the extra dimensions remained 
comparatively small during cosmological evolution. 

Such a mechanism was proposed in [24, 25] . The basic 
idea is that strings dominate the dynamics of the early 
universe and can see each other most efficiently in 2(p+l) 
(p=l for strings) dimensions. Therefore, strings can only 
interact in three spatial dimensions, while strings mov- 
ing in higher dimensions eventually cease to interact effi- 
ciently and their winding modes will prevent them from 
further expanding. If branes are included, it was shown 
in [26] that strings will still dominate the evolution of 
the universe at late times so the mechanism of [24] still 
survives. 



In six dimensions and for codimension-2 braneworlds, 
the gravity dynamics appear even more radical and still 
a good understanding of cosmology and more generally 
gravity in such theories is missing. The most attractive 
feature of codimcnsion-2 braneworlds is that the vacuum 
energy (tension) of the brane instead of curving the brane 
world-volume, merely induces a deficit angle in the bulk 
solution around the brane [27] (see also [28] for string-like 
defects in six dimensions). This looks very promising in 
relation to the cosmological constant problem, although 
in existing models, nearby curved solutions cannot be ex- 
cluded [29] , unless one allows for singularities more severe 
than conical in particular supersymmetric models [30] . It 
was soon realized [31] that one can only find nonsingular 
solutions if the brane energy momentum tensor is pro- 
portional to its induced metric, which means simply that 
it is pure tension. A non-trivial energy momentum ten- 
sor on the brane causes singularities in the metric around 
the braneworld which necessitates the introduction of a 
cut-off (brane thickness) [32-34]. 

An alternative approach to study the gravitational dy- 
namics of matter on infinitely thin branes is to modify 
the gravitational action as discussed previously. Indeed, 
it was shown in [35] that the inclusion of a Gauss-Bonnet 
term in the gravitational action allows a non-trivial en- 
ergy momentum tensor on the brane, and in the thin 
brane limit, four dimensional gravity is recovered as the 
dynamics of the induced metric on the brane. The pe- 
culiar characteristic of this way to obtain four dimen- 
sional gravity for codimension-2 branes, is that, apart 
from the inclusion of a (deficit angle independent) cos- 
mological constant term, there appear to be no correc- 
tions to the Einstein equations coming from the extra 
dimensions in the purely conical case. Another possi- 
bility, discussed in [36], is to study (instead of conical 
3-branes) codimension-2 branes sitting at the intersec- 
tion of codimcnsion-1 branes in the presence again of a 
bulk Gauss-Bonnet term. 

Much less has been done, however, for the cosmology of 
theories in six or higher dimensions with branes of codi- 
mension greater than one. This is because, unlike the 
codimension one case, these branes exhibit bulk curva- 
ture singularities which are worse than (5-function singu- 
larities. They then need some regularization (introduc- 
tion of brane thickness) which makes the study of cos- 
mology on them rather complicated [37]. An alternative 
way to study cosmologies of branes of higher codimen- 
sion would be to consider corrections to the gravitational 
action, such as an induced curvature term on the brane 
[14] and a Gauss-Bonnet term in the bulk [17], which al- 
low the brane to have a mild singularity structure (see 
also [39]). These thin brane cosmologies would have the 
additional advantage that the internal structure of the 
brane does not influence the macroscopic cosmological 
evolution. 

Six-dimensional models can also be obtained in super- 
gravity theories. In six-dimensional supergravity theory 
there is a vacuum solution which has the structure of 
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Minkowskix5 2 [40, 41]. Cosmological applications of 
this solution was presented in [42]. These theories are 
interesting because many of the parameters of the mod- 
els are fixed from the theory and in particular the form 
of the scalar potential is dictated from the structure of 
the supergravity theory. The drawbacks of these models 
however is that they are plagued from anomalies. The 
reason is that they are chiral and therefore a gauging 
of some symmetries is required for cancellation of these 
anomalies [41]. 

In this talk we will discuss braneworld models of 
codimcnsion-1 and codimension-2 in six dimensions ana- 
lyzing their cosmology. We will also briefly discuss cos- 
mological models in six-dimensional supergravity. First 
we will discuss a (4 + 1) braneworld model in a six- 
dimensional spacetime bulk with a cosmological constant. 
We will consider a 4-brane fixed at some position in 
the sixth dimension and we derive the dynamical six- 
dimensional bulk equations in normal gaussian coordi- 
nates with a cosmological constant in the bulk and con- 
sidering matter on the brane [43, 44]. We look for time 
dependent solutions allowing for two scale factors, the 
usual scale factor a{t) of the three dimensional space 
and a scale factor b(t) for the extra fourth dimension. 
If a(t) — b(t) we get the Friedmann equation of the gen- 
eralized Randall- Sundrum model in six dimensions de- 
scribing a four-dimensional universe. If a(t) ^ b(t) we 
get a generalized Friedmann equation in six dimensions. 

The problem can be looked at a different angle [45]. 
If a = b, we can write the bulk six-dimensional met- 
ric in "Schwarzschild" coordinates and then we have the 
equivalent description of a 4-brane moving in a static 
Schwarzschild- (A) dS six-dimensional bulk. If however 
a 7^ b, a brane observer uses a and b, for whom they 
are static quantities, to measure the departure from six- 
dimensional spherical symmetry of the bulk. The im- 
portant result of this consideration is that [46] , since the 
brane observer needs to define a cosmic time in order to 
derive an effective Friedmann equation on the brane, a 
and b are related through the Darmois-Israel junctions 
conditions and because of that, their relation depends 
on the energy-matter content of the brane. The physical 
reason of the existence of such a relation between a and 
b is that the requirement of having a cosmological evolu- 
tion on the brane introduces a kind of compactification 
on it and the relation between a and b acts as a constraint 
of the brane motion in the bulk. 

Using the six-dimensional generalized Friedmann equa- 
tion and assuming that p ^ p, where p is the pressure 
of the physical three dimensions and p corresponds to 
the fourth dimension, we make a systematic numerical 
study of the cosmological evolution of the scale factors 
a(t) and b(t) for several values of the parameters of the 
model, A@ the six-dimensional cosmological constant, k 
the brane spatial curvature and w and w parameteriz- 
ing the form of the brane energy-matter content of the 
three dimensions and of the extra fourth dimension re- 
spectively. We find that [46] , in order the fourth dimen- 



sion to be small relatively to the other three dimensions 
and to remain constant during cosmological evolution, w 
must be negative, indicating the presence of a dark form 
of energy in the extra fourth dimension. We find this 
result for all cases considered, (A)dS or Minkowski bulk, 
open, closed or flat universe, radiation, dust, cosmologi- 
cal constant and dark energy dominated universe and the 
specific value of to depends on the energy-matter content 
of the other three dimensions. 

A codimension-2 braneworld model in six dimensions 
will be discussed next. Results similar to the Gauss- 
Bonnet case [35], i.e. four dimensional gravity for an ar- 
bitrary energy momentum tensor, can be obtained if we 
include in the action an induced gravity correction term 
instead. Again, in the purely conical case, there appear 
to be no corrections to the Einstein equations coming 
from the extra dimensions. The most important obser- 
vation is that the brane and bulk energy momentum ten- 
sors are strongly related and any cosmological evolution 
on the brane is dictated by the bulk content [48] . We also 
see how this correlation is relaxed in the case where bulk 
Gauss-Bonnet terms and brane induced gravity terms are 
combined. Thus, the necessary presence of extra curva- 
ture terms in the gravitational action in order to give non- 
trivial gravitational dynamics on a codimension-2 brane, 
leads to a realistic cosmological evolution on the brane in 
the thin brane limit, only if a Gauss-Bonnet term is in- 
cluded. However, let us note that the most physical way 
to investigate the dynamics of codimcnsion-2 branes, is 
by giving thickness to the brane [33, 34]. 

We will discuss in details the cosmological evolution 
of a conical brane with both an induced gravity and 
a Gauss-Bonnet term added in the higher dimensional 
gravity action [49]. For simplicity we will assume that 
the only matter in the bulk is cosmological constant A b ■ 
We will solve the equations of motion evaluated on the 
brane and assume that the integration of them in the 
bulk does not give rise to pathologies (i.e. singularities). 

Firstly, we will discuss the isotropic cosmology, in 
which the brane matter has to obey a tuning relation. 
The physical meaning of this relation is that for any mat- 
ter wc put on the brane its " image" should be present in 
the bulk. We will see that the evolution of the system 
for Kb = tends to a fixed point with to = 1/3 and for 
Ab > to a fixed point with w = —1. For As < the 
system has a runaway behaviour to w — > oo. 

We will then relax the isotropy requirement for 
the metric (keeping, however, the matter distribution 
isotropic) in order to find whether the above matter tun- 
ing is an attractor or not. The matter on the brane need 
not now satisfy the previous tuning relation and the al- 
lowed regions of initial values of the energy density and 
pressure are significantly larger. The analysis of the dy- 
namics of the system shows that line of isotropic tuning is 
an attractor for A^ > and thus the system isotropises 
towards it. However, for values of A^ which give a real- 
istic cosmological evolution of the equation of state, the 
attractor property of the previous line is very weak and 
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fine tuning of the initial conditions is necessary in order 
to have an evolution with acceptably small anisotropy. 
For Ab < the system shows, as in the isotropic case, a 
runaway behaviour. 

Finally we will discuss in brief, six-dimensional cosmo- 
logical models coming from supergravity theories. These 
models arc interesting because, they give conventional 
four-dimensional cosmology with less arbitrariness in cos- 
mological parameters and less fine-tuning. 



2. A (4 + 1) COSMOLOGICAL 
BRANEWORLD MODEL IN 
SIX-DIMENSIONAL SPACETIME 

In this section we describe a braneworld cosmological 
model of codimension-1 in six dimensions. This model 
is a generalization of the existing braneworld models in 
five-dimensions. We study the model following two dif- 
ferent approaches. First we use normal gaussian coor- 
dinates to describe a static 4-brane at fixed position in 
a six-dimensional spacetime bulk. The six-dimensional 
Einstein equations are derived and with the use of the ap- 
propriate junction conditions the generalized Friedmann 
equation is given. The cosmological evolution described 
by this generalized Friedmann equation involves the usual 
three-dimensional scale factor a(t) and the scale factor 
b(t) describing the cosmological evolution of the extra 
fourth dimension. 

We consider next a dynamical brane moving in a bulk 
described by six-dimensional static "coordinates" . In this 
case the dynamical brane is moving on a geodesic which 
is given by the junctions conditions. We derive the equa- 
tions of motion of the brane which for a brane observer 
describe the cosmological evolution on the brane. We 
discuss the connection between the static and dynamical 
brane models and the physical information that can be 
extracted from these approaches. 



2-1. Static Brane in a Dynamical Bulk 

We look for a solution to the Einstein equations in six- 
dimensional spacetime with a metric of the form 

ds 2 = -n 2 (t,y 1 z)dt 2 + a 2 (t 1 y,z)dYj 2 k 

+ b 2 (t,y,z)dy 2 + d 2 (t,y,z)dz 2 , (2.1) 

where dT, 2 . represents the 3-dimcnsional spatial sections 
metric with k = —1, 0, 1 corresponding to the hyper- 
bolic, fiat and elliptic spaces, respectively. 

If the brane is fixed at the position zo, then the to- 
tal energy-momentum tensor can be decomposed in two 
parts corresponding to the bulk and the brane as 



where the energy-momentum tensor on the brane is 

T M(6) = fcp) diag{ _ pMp )iJj0)) (2 . 3) 

where p is the pressure in the extra brane dimension. 
We assume that there is no matter in the bulk and the 
energy-momentum tensor of the bulk is proportional to 
the six-dimensional cosmological constant. 

The presence of the brane in z imposes boundary con- 
ditions on the metric: it must be continuous through 
the brane, while its derivatives with respect to z can be 
discontinuous at the brane position. This means that 
the generated Dirac delta function in the metric sec- 
ond derivatives with respect to z must be matched with 
the energy-momentum tensor components (2.3) to satisfy 
the Einstein equations. The Darmois-Israel conditions 
are [46], 
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(p - p + p) , 



{ P -3(p-p)} , 



(2.4) 
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{p + 3{p + p)} , 



where the subscript (0) indicates quantities on the brane. 
The energy conservation equation on the brane can be de- 
rived taking the jump of the (06) component of the Ein- 
stein equations and using the junction conditions (2.4) 
and the corresponding time derivatives. We obtain 



p + 3(p + p)^ + (p + p)^ 



0. 



(2.5) 



To find the Friedmann equation we take the jump of 
the (66) component of the Einstein equations and we use 
the fact that 



[d z fd z g] = #d z f# [d z g] + [d z f] #d z 9# ■ 



where 



#/(«)# : 
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is the mean value of the function / through y 
we arrive to the following equation 



-v = r - 



3 P - 



(2.6) 



0, and 



(2.7) 
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Taking the mean value of the (66) component of the 
Einstein equations and using the junction equations, 
(2.7) (we have assumed a Z 2 symmetry), we arrive to 
the generalized Friedmann equation in six-dimensions 

/ ao 1 bo apbo a| \ _ 
\a 3b a b Q j 



32 



p(p + 2p+ |p) + (p - pf 



rpM _ rfiM(B) rpM(b) 

1 N — 1 N N ) 



(2.2) 



a 2 U 2 ^' 



(2.8) 
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where we have assume (d y a) = 0, (d 2 a) = and we 
have chosen n = 1. 

In the case of a 3-brane in a five-dimensional bulk, the 
first integral of the space-space component of the Ein- 
stein equations, with the help of the other equations can 
be done analytically [9] and this results in the Friedmann 
equation on the brane with the dark radiation term as an 
integration constant [50]. In our case the Einstein equa- 
tions cannot be integrated analytically and therefore, the 
usual form of the Friedmann equation on the brane can- 
not be extracted from (2.8). Nevertheless, if a(t) and 
b(t) are related, this equation will give the cosmological 
evolution of the scale factor a(t). 

If a(t) = b(t) = 1Z(t) then (2.8) becomes 
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(2.9) 

This equation is the generalization of the Randall- 
Sundrum Fricdmann-likc equation in six dimensions, and 
it has, as expected, the p 2 term with a coefficient adjusted 
to six dimensions. Note that this equation can easily be 
generalized to D dimensions. 



2-2. Dynamical Brane in a Static Bulk 

We consider a 4-brane moving in a six-dimensional 
Schwarzschild-AdS spacetime. The metric in the 
"Schwarzschild" coordinates can be written as 

z 2 

ds 2 = -h(z)dt 2 + —dY? k + h- 1 (z)dz 2 , (2.10) 



where 



dZ 2 = 



dr 2 



1 — kr 2 



r 2 dVL 2 {2) + (1 - kr 2 )dy 2 , (2.11) 



and 



h{z) = k 



I 2 



(2.12) 



Comparing the metric (2.1) with the metric (2.10) we can 
make the following identifications 



n(z) 
a(z) 
d(z) 



b(z) = z/l 



(2.13) 



The two approaches of a static brane in a dynamical bulk 
described in Sec. 2. 2-1 and of a moving 4-brane in a static 
bulk are equivalent [51]. To prove this equivalence, con- 
sider the Darmois-Israel conditions for a moving brane, 
[K£\ = -k 2 6) (T>? - \Th$), with h*t being the induced 
metric on the brane; analogously to the case in 5 dimen- 
sions one can obtain 



-ti - 2TZ 
Vh + TZ 2 



(2.14) 
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Vh + TZ 2 «( 6) 



K 



P, 



(2.15) 



where we have defined z = TZ(t). Combining these two 
equations and using (2.12) we find that 

2 | + 3 f|V = -3^p 2 -^-3A_A. (2 . 16 ) 
1Z \K 64 1 8 1Z 2 I 2 ' 



Comparing with (2.9) we see that 
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(2.17) 



with Tg = — A 6 the bulk cosmological constant, and I the 
size of the AdS space. Therefore, a brane observer de- 
scribes the cosmological evolution of a four-dimensional 
universe with the Friedmann equation (2.16) with 1Z pa- 
rameterizing the motion of the brane in the z direction. 

To describe the motion of the 4-brane with a(z) ^ b(z) 
we can write the metric (2.10) as 

ds 2 = -n(z) 2 dt 2 + a 2 {z)dY,l + b(z) 2 dy 2 + d 2 (z)dz 2 . 

(2.18) 

We denote the position of the brane at any bulk time t 
by z = TZ(t) as before. Then, an observer on the brane 
defines the proper time from the relation 



n 2 (t,1Z{t))t 2 - d 2 {t,1Z(t))iZ 2 = 1 , 



(2.19) 



which ensures that the induced metric on the brane will 
be in FRW form 



ds ir 



dr 2 

,2 



- n 2 (t,1Z(t))i 2 - d 2 (t,TZ(t))lZ 2 
+ a 2 {t,1Z{t))dY? (3) +b 2 (t,1Z(t))dy 2 
-dr 2 +a 2 (t,1Z{t))dY l 2 3) 
b 2 {t,1Z{t))dy 2 , (2.20) 



where the dot indicates derivative with respect to the 
brane time r. 

Introducing an energy-momentum tensor on the brane 



J- flU ll VCL ± fj, 



1, 



(2.21) 



where = diag(—p,p,p,p,p), the Darmois-Israel con- 
ditions become 



" K (6)^V" 



(2.22) 



where K^ v is the extrinsic curvature tensor. These give 
the equations of motion of the brane 



d 2 dV? - dlZ yj + d 2 TZ 2 
Vl + d 2 1Z 2 n 



(dnK+^ 
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d z a 
ad 

dzb 
bd 



y/l + d?li 2 
Vl + d 2 ^ 



— (dd z n — nd z d)1Z 2 ^ 

= -^(3(p + p)+p) , (2.23) 

= -t(' + '-')' 
= j|>(p-3(p-fl) 



(2.24) 
(2.25) 



Notice that if a = 6 = we recover equations 

(2.14) and (2.15), corresponding to the Schwarzschild- 
AdS spacetime. 

Equation (2.23) is the main dynamical equation de- 
scribing the movement of the brane-universe in the six- 
dimensional bulk, while a combination of (2.24) and 
(2.25) acts as a constraint relating a and b (remember 
that for a brane observer a and b are static, depending 
only on z) 



j{ b(p+P-p)/(p-3(p-p)) 



(2.26) 



where A is an integration constant. This means that the 
relative cosmological evolution of a, the scale factor of 
the three-dimensional physical universe and b, the scale 
factor of the extra dimension, depends on the dynamics 
of the energy-matter content of the brane-universe. 

Another interesting observation is that, a brane 
observer can measure the departure from full six- 
dimensional spherical symmetry of the bulk using the 
quantities a and b. If a(z) = b(z) = z/l then the 
symmetry of the bulk is S 4 having a six-dimensional 
Schwarzschild-(A)dS black hole solution. If a{z) ^ b(z) 
fixing a to be a(z) — z/l, because of (2.26) b(z) is given 

by 



b= 4 



l_,z 
A\l 



(p-3(p-p)) / '(p+p-p) 



(2.27) 



and we expect the topology of the bulk to be S 3 x A4, 
M being a compact or non-compact manifold. There are 
no analytical solutions in the six-dimensional spacetime 
with such topology and the reflection of this on the brane 
is the difficulty of the brane equations (2.23)-(2.25) to be 
integrated analytically. Note also that because of (2.27), 
a change in the topology of the bulk is triggered by the 
dynamics of the matter distribution on the brane. 

It is also interesting to further explore the rela- 
tion we found between the cosmological evolution of a 
higher-dimensional brane-universe with the static prop- 
erties of the bulk. The cosmological evolution on a 
higher-dimensional brane-universe is related to a topol- 
ogy change of the bulk during the evolution, and this rela- 
tion might lead to a better understanding of the Gregory- 
Laflamme [47] instabilities of higher-dimensional objects. 



3. THE COSMOLOGICAL EVOLUTION 
OF A (4+1) BRANE-UNIVERSE 

To study the cosmological evolution of the four- 
dimensional brane-universe, we made the following as- 
sumptions for the initial conditions and the matter dis- 
tribution on the brane. We assume that the universe 
started as a four-dimensional one at the Planck scale, all 
the dimensions were of the Planck length and the matter 
was isotropically distributed. In this case the cosmolog- 
ical evolution is described by the generalized Friedmann 
equation (2.9). Then an anisotropy was developed in the 
sense that p = Qp with Q ^ 1. The cosmological evo- 
lution is now described by (2.8) supplemented with the 
constrained equation (2.26) and the matter distribution 
on the brane is given by the equations of state 



p = wp 
p = wp 



(3.1) 
(3-2) 



Using (2.26) and relations (3.1) and (3.2) the generalized 
Friedmann equation (2.8) becomes 

w-w ■ (1 + w)(-3w + 2w - 1) + 3(1 + w) 2 2 

H a + — H„ 



1 + w 



+ 



(1+w) 2 
2w - w + 1 p „ 2 
(1 + w) 2 ~p 
P 



H„ 



3(1 + w) 2 p 2 



3(l + M>)/9 



= -^-\l + 2w+-w + (w-wf 



+ 



"(6) 



A fi . 



(3.3) 



Using the conservation equation (2.5) to eliminate p and 
its derivatives from the above equation, the cosmologi- 
cal evolution of the three-dimensional scale factor a(t) is 
given by the equation 



da 2C+1 + 



+ 



B 2 2B 
32 



l + 2w + -w + (w 



— a 



2C 



a 2 ^A 6 -2fc 



= 



where the constants B and C are given by 
1 - 3w + 3w 

B 

c 



1 + w — w 
3(1 + w) +B(w + 1) 



while the b(t) scale factor is 

b(t) = a(t) B 



wf 
(3.4) 

(3.5) 
(3-6) 

(3.7) 



We made a numerical analysis of equations (3.4) and (3.7) 
and studied the time evolution of the two scale factors for 
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different backgrounds and spatial brane-curvature. We 
allowed for all possible forms of energy-matter on the 
physical three dimensions (w=0, 1/3, -1/3) and also for 
the possibility of dark energy (w=-l), leaving w as a free 
parameter. 

For the scale factor b(t) to be small compared to the 
scale factor a(t), the constant B in (3.7) should be nega- 
tive. In Table 1 we give the allowed range of values of w 
for various values of w. These values in turn were used to 
plot the time evolution of the scale factors a(t) and b(t) 
using (3.4) and (3.7) respectively. The criterion for the 
acceptance of a solution is to give a growing evolution of 
a(t) and a decaying and freezing out evolution for b(t). 
The results for various choices of the parameters of the 
model are presented in the following table 



w 


w 


-1 


> or < -4/3 





> 1 or < -1/3 


1/3 


> 4/3 or < 


-1/3 


> 2/3 or < -2/3 



TABLE I: The allowed values of w for B to be neg- 
ative. 

The numerical analysis showed [46] that in all the cases 
considered, w is negative in the range of values given in 
Table 1 for all acceptable solutions, indicating the need 
of dark energy to suppress the extra fourth dimension 
compared to the three other dimensions. 



4. CODIMENSION-2 BRANEWORLD 
MODEL WITH INDUCED GRAVITY IN 
SIX DIMENSIONS 

We will next discuss a codimension-2 braneworld 
model in six dimensions. We will first consider a six- 
dimensional theory with general bulk dynamics encoded 
in a Lagrangian Lsuik and a 3-brane at some point r = 
of the two-dimensional internal space with general dy- 
namics Lbrane in its world-volume. If we include an 
induced curvature term localized at the position of the 
branc, the total action is written as: 



Mi 



d 6 xVGR {6) + . 



***** 



+ 



J d 6 xC B uik + J d*x£ 



brane 



J(r) 
2ttL 



(4.1) 



In the above action, M 6 is the six-dimensional Planck 
mass, M 4 is the four-dimensional one and r c — M^/Mq 
the cross over scale between four-dimensional and six- 
dimensional gravity. The above induced term has been 



written in the particular coordinate system in which the 
metric is 



dsl 



g^ u (x,r)dx^dx v + dr 2 + L 2 (x,r)d9 2 



(4.2) 



where g tlv {x 1 0) is the braneworld metric and x^ denote 
four non-compact dimensions, jj, = 0,...,3, whereas r,9 
denote the radial and angular coordinates of the two ex- 
tra dimensions (the r direction may or may not be com- 
pact and the 9 coordinate ranges form to 2tt). Capi- 
tal M,N indices will take values in the six-dimensional 
space. Note, that we have assumed that there exists an 
azimuthal symmetry in the system, so that both the in- 
duced four-dimensional metric and the function L do not 
depend on 9. The normalization of the <5-function is the 
one discussed in [52]. 

To obtain the braneworld equations we expand the 
metric around the brane as 



L{x,r) = (3{x)r + 0{r 2 ) 



(4.3) 



At the boundary of the internal two-dimensional space 
where the 3-brane is situated, the function L behaves as 
L'(x,0) — (3(x), where a prime denotes derivative with 
respect to r. As we will see in the following, the de- 
mand that the space in the vicinity of the conical singu- 
larity is regular, imposes the supplementary conditions 
that dft/3 = and d r g l _ tv (x, 0) = 0. 

The Einstein equations which are derived from the 
above action in the presence of the 3-brane are 



M " v 2irL 



Mi 



T (B)N T (br)v?n rN ( H r ) 

Hi +H "m", 2nL 



(4.4) 



with G { $ N and Gp )v the six-dimensional and the four- 
dimensional Einstein tensors respectively, the bulk 
energy momentum tensor and y^ fcr ) 1/ the brane one. 

We will now use the fact that the second derivatives 
of the metric functions contain (5-function singularities at 
the position of the brane. The nature of the singularity 
then gives the following relations [35] 

L" S(r) 

—— = —(1 — L')— h non — singular terms (4.5) 

K' v S(r) 

—f— = K^ v — h non — singular terms . (4.6) 

From the above singularity expressions we can match the 
singular parts of the Einstein equations (4.4) and get the 
following "boundary" Einstein equations 

1 9ir 9ir T 



r 2 M 6 



(4-7) 

where K^ u is the extrinsic curvature and we denote by 
|o the value of the corresponding function at r — 0. 
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We will now make the assumption that the singularity 
is purely conical. In the opposite case, there would be 
curvature singularities R^ oc 1/r, because in the Ricci 
tensor R,ffi there are terms of the form [35] 



»(«) - --—d a 



a, 



2r 



0(1) , (4.8) 



which in the vicinity of r = arc singular if d r g^ v {x^ 0) ^ 
0. The absence of this type of singularities imposes the 
requirement that K^\o = 0. Then, the Einstein equa- 
tions (4.7) reduce to 



1 9tt 

T^ + ^(l-[3) 9 ^\ . (4.9) 



r*M$ 



The four-dimensional Einstein equations (4.9) describe 
the gravitational dynamics on the brane [48]. The effec- 
tive four-dimensional Planck mass and cosmological con- 
stant arc simply 



M 2 Pl = M 2 A = r 2 M? 



4 - < c -6 . (4-10) 
A 4 = A - 2ttM 6 4 (1 - p) , (4.11) 



where A is the contribution of the vacuum energy of the 
brane fields. The normalization of A4 is defined by the 
convention that the four dimensional Einstein equation 

reads Guv = -tj^-{T U u — A4gn U ). Note that in contrast 

m pi 

to the case of [35], the four dimensional Planck mass is 
independent of the deficit angle. 

Furthermore, it is interesting to see that contrary to 
the five-dimensional case, the induced gravity term in 
six dimensions does not introduce any correction terms, 
apart from a cosmological term, in the four-dimensional 
Einstein equations on the brane, unless singularities of 
other type than conical are allowed in the theory, and 
a regularization scheme is employed. In the latter case, 
the last term of the right hand side of (4.7) would provide 
information of the bulk physics. This absence of correc- 
tions in the purely conical case, is exactly what happens 
also in the case of the bulk Gauss-Bonnet theory [35] . 

What is important to note at this point, is that al- 
though we have found a "boundary" Einstein equation, 
there is more information about the dynamics of the the- 
ory contained in the full six-dimensional Einstein equa- 
tions. In [48] it was showed that the energy momentum 
tensor of the bulk is strongly related to the energy mo- 
mentum tensor of the brane via the relation 



1 

27? 



-8ttM 6 4 (1-/3) 



(4.12) 



This equation constitutes a very strong tuning relation 
between brane (xf r)A1 ) and bulk (ri B)r |o) matter. It 
shows that, in order to have some cosmological evolu- 
tion on the brane (i.e. time dependent T^) and since (5 
is constant, the bulk content should evolve as well in a 
precisely tuned way. 



We can compare the above result with what is hap- 
pening in five dimensions. The action of a general five- 
dimensional theory with a 3-brane at the point r = of 
the extra dimension, and with an induced curvature term 
localized on it, 



is 



o _ Mi 



J dPxVGRW +r c j d 4 x^gR {4) £(: 
d 5 x£ B uik + / d 4 x£ brane S(r) . 



(4.13) 



In the above action, M5 is the five-dimensional Planck 
mass, M4 is the four dimensional one and r c = M 2 /M^ 
the cross over scale of the five-dimensional theory. The 
above induced term has been written in the particular 
coordinate system in which the metric is written as 



ds 5 = g ul/ (x,r)dx ,J 'dx v + dr 2 



(4.14) 



where the x M denote the usual four non-compact dimen- 
sions, fi = 0, ...,3, whereas r denotes the radial extra 
coordinate and capital M,N indices will now take values 
in the five-dimensional space. 

The Einstein equations which are derived from the ac- 
tion (4.13) in the presence of the 3-brane are 



Ml 



T^ N + T^>6^S(r) , (4.15) 



with the bulk energy momentum tensor and 

the brane one. Performing a similar analysis like the six- 
dimensional case [48] we get 



r(% + \(k;kz 



K 2 



Ml 



T^ r \ . (4.16) 



From the above equation we see that the bulk matter 
content does not necessarily dictate the brane cosmolog- 
ical evolution. This is because the extrinsic curvature on 
the brane can be non-trivial and it is this one which 
plays the most crucial role in the cosmology. In other 
words, in five dimensions it is the freedom of the brane 
to bend in the extra dimension which makes the evolu- 
tion not tuned to the bulk matter content. The absence 
of such bending in six dimensions (imposed by singular- 
ity arguments) gives the bulk the crucial role for how the 
brane evolves. 

We can also introduce a Gauss-Bonnet term in the six- 
dimensional action (4.1) and see how the above results 
are modified. In this case the action (4.1) is augmented 
by the term 



Sgb 



M%a 



c^ (6)2 -4<y + <y KA ). (4.17) 



Then the variation of the above action introduces an ex- 
tra term in the left hand side of the Einstein equations 



-9- 



(4-4), 



H 



AR 



2R^R% )N 



(6) 2 
KA 

NP 



d(6) 2 x 

n-ABKk) 



, /| p( 6 ) N -T.KP op(6) r>NKAP 
+ W KMP U (6) ~ ZU MKAP U (6) 



. (4.18) 



Equating the singular terms of the Einstein equations by 
the standard procedure of section 2, and demanding that 
the singularity is purely conical, we obtain the following 
"boundary" Einstein equations 



1 



+ 



M4(r 2 + 8tt(1 
2tt(1 - g) 
r 2 + 8tt(1 - /3)a 



T (br) 

d)a) 



(4.19) 



Equation (4.19) describes the gravitational dynamics on 
a codimension-2 brane when both induced gravity and 
Gauss-Bonnet correction terms are present. The effective 
four-dimensional Planck mass and cosmological constant 
are simply 



M 



,., -= Ml(r 2 c + 8tt(1 - p)a) , (4.20) 
= A-2ttM 6 4 (1-/3) , (4.21) 



where A is the brane tension. Note that the Planck mass 
this time can depend on the deficit angle. This is an 
effect of solely the bulk Gauss-Bonnet term. 

Evaluating the the (rr) component of the Einstein 
equation at the position of the brane r = we obtain 
the following relation [48] 



^ (4) lo + a(R^ 2 - AR% 2 + <L 2 A )| = -^T<*> 



o • 



(4.22) 

From (4.22) we see that there can be no relation be- 
tween the extra dimensional component Tr B ^ r \o of the 
bulk energy momentum tensor at the position of the 
brane with the brane energy momentum tensor X^ 6 ^. 
This is due to the appearance of the Riemann curvature 
which cannot be evaluated from previous equations (the 
Ricci tensor and scalar can be substituted from (4.19)). 

Instead, using (4.19) and (4.22), one can solve for R^^ x \o 
as a function of the brane and bulk matter at the position 
of the brane. 



5. COSMOLOGICAL EVOLUTION OF A 
CONICAL CODIMENSION-2 
BRANEWORLD MODEL 

We consider a six-dimensional theory with general bulk 
dynamics encoded in a Lagrangian Csuik and a 3-brane 
at some point r = of the two-dimensional internal space 
with general dynamics Lbrane in its world-volume. The 



gravitational dynamics is described as we discuss in the 
previous sections by a Gauss-Bonnet term in the bulk and 
an induced four-dimensional curvature term localized at 
the position of the brane. Then the total action is written 
as 



c _ Mi 



- AR 



+ 



{ J d e xy/^)[R^ + a(R ( V 2 
i} 2 + R {6) 2 )] + r 2 [ tfxsf^R^ 

IN + n MNKA>\ + r c / « X V 9 n 2irL) 

/f 6(r) 
d 6 x£ Bulk + / d 4 xC brane ^—j- . (5.1) 



The full equations of motion that are derived from the 
above action using the metric (4.2) are 



M 6 ) N _1_ r 2 r (4WM xN^j)_ ttN 



1 



T (B)N T (br)vw ,N 6 ( r ) 
± M ' 1 fJ. °M°v 2 n £ 



(5.2) 



with G 



(6) N 
M 



R 



(6) N 
M 



Einstein tensor, G^ v = R^ v 



\R^&m the six-dimensional 



\R {A) 8^ the four- 



dimensional Einstein tensor and is given by (4.18). 

In order that there are no curvature singularities more 
severe than conical, we will impose certain conditions on 
the value of the extrinsic curvature = g' on the 
brane, where the prime denotes derivative with respect 
to r, and on the expansion coefficients of the function L 



(3i(x)r + /3 2 (x)r 2 + f3 3 (x)r 3 



These conditions read [35] 
[3\ = const. 



= 0, 

and p2 = 



(5.3) 



(5.4) 
(5.5) 



Imposing these conditions and keeping only the finite 
part in L"/L, the Einstein equations (5.2) can be eval- 
uated at r — 0. The effective Einstein equations on the 
brane (obtained by equating the (5-function parts of the 
Einstein equations) are 



" M 2 



pi 



T^ r > - A 4 <^ 



(5.6) 



with M 2 pl = M 4 [r 2 c + 8tt(1 - ft) a] and A 4 = -2ttM|(1 - 
ft ) . The various components of the bulk Einstein equa- 
tions evaluated at r — are given in the following: 
The (fj,v) component 



—a 



+AK? R> 



L" 

2K'R-A — R 
L 



—2RR" + AR UK R U 
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9 7? ryVKXp 

+4^-r; + 2k'r; + k;'r 



-2K: 'R; 2g, x K^ 'R™\ = ^Tf K (5.7) 



The (rr) component 
2 2 V 



_J_ T (B)r 

Mr r 



(5i 



The (09) — (rr) component 
\K' + a (K'R 2K? 'R;) = -^(Tl B)e t£ 
The (fi, r) component 



B)r\ 



T^ r = . 



(5.9) 



(5.10) 



In the following, we will study the above equations in a 
time dependent background, assuming that the bulk con- 



sists of a pure cosmological constant T 



(B)N 
M 



B°M 



-Ab5 



and that the matter content of the brane is an isotropic 
fluid with 1$ r)v = 5% diag(-p 6 , P b , P b , P b ). 

5-1. The Constrained Isotropic Case 

We are interested in the cosmological evolution of a 
flat isotropic brane-universe, therefore we will consider 
the following time dependent form of the metric (4.2) 

ds 2 = -N 2 (t,r)dt 2 + A 2 (t,r)dx 2 + dr 2 + L 2 (x,r)d9 2 . 

(5.11) 

We can use the gauge freedom to fix N(t,0) = 1, while 
we define A(t, 0) = a(t). The curvature singularity avoid- 
ance condition (5.4) we imposed, dictates that N'(t, 0) = 
A'(t, 0) = 0, while the second derivatives of these metric 
functions are unconstrained. 

For this ansatz, the Einstein equations (5.6) and (5.8) 
give 

(5.12) 
(5.13) 

w< ■ < 5 -»> 

The equations (5.12) and (5.13) are the usual Fried- 
mann and Raychaudhuri equations of a four-dimensional 
universe with a scale factor a, while the third equation 
(5.14) appears because of the presence of the bulk and 
acts as a constraint between the matter density and pres- 
sure on the branc. To see this, a simple manipulation of 





4 

er 


Pb + A 4 

M 2 
ivi pl 


2*- 


a 2 


-P b + A 4 


a 


f — 


Mh 


d 


d 2 \ 


aa 2 


- + 




+ 12a— = 


a 


a 2 


or 



the above equations gives 



Mi 



- 2 



1 
2 " 

A 4 
M 2 PI 



2 A 4 



pi 



iPb - Pb 



. 2 A 4 
1+ 3 a Ml 



pi 



2 p b (3P b - 
3 a M 2 



Pb) 



(5.15) 



pi 



which shows a precise relation between As, p b and P b . 
To simplify the equations, we assume that the vacuum 
energy (tension) of the brane cancels the contribution A 4 
induced by the deficit angle, i.e. 

p b = -A 4 + p m and P b = A 4 + P m , (5.16) 

with P m = w c p m . Then the above equations read 

Pm 



4 



'Ml, 



(5.17) 
(5.18) 



while the constraint equation becomes 
A B p 



Ml 



i) + |(3«, c + i)a£| 



pi 



(5.19) 

Using the metric (5.11), the Einstein equations (5.7) 
and (5.9) can be solved for the second r-derivatives of the 
metric, as functions of the matter content on the branc, 
and in principle they can give us information about the 
structure of the bulk at r = 



A^_ 
~A 

~N~ 
T 



l(l + 4aA|) K + l)^ 



Mr. 



= 3- 



A I 



w c - 1 

1 



(w c + l)(3w c -l) 



4a 



1 



1 + 4q;-£is 
~A 



As 

Ml 



A" / 

-3 — \l + 2a(w c + l) 



(5.20) 



A potential problem in the cosmology of the system 
would be, if the denominator of (5.20) is equal to zero, 



i.e. when w c - 



We with 



4a 



± 



64a^ 



32a^ 



13 



12a^ 



(5-21) 

When this happens, the six-dimensional curvature invari- 
ant will diverge close to the brane. Thus, after discussing 
the cosmological evolution of the brane world- volume, we 
should always check that it does not pass through a point 
which satisfies the above relation. 
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5-2. Cosmological Evolution of an Isotropic 
Brane-Universe 

From the constraint relation (5.19), we can solve for 
w c , the allowed equation of state of the matter on the 
branc. It should satisfy the following equation 



-2 



+ 



Pm 

Ml 



(5.22) 



with p m > so that the Hubble parameter is real. 

Before analyzing the system, let us note a first im- 
portant difference between the system of the pure four 
dimensional dynamics and the one with the extra con- 
straint added because of the presence of the bulk. In 
the four dimensional system, a constant w is allowed and 
its value is preserved during the evolution of the system. 
On the other hand, the evolution of the system with the 
extra constraint forbids any evolution with constant w c . 
Indeed, by differentiating (5.22) and using (5.17) and the 
conservation equation p m + 3(1 + w c )p m ^ = 0, we can 
find a differential equation for w c 



W c + 3(1 + W c )p r 



dw c a 
dp m a 



. 



(5.23) 



Imposing a further condition of keeping w c constant, 
would result to a constant p m , related to w c in a specific 
way, and by the conservation equation, to zero Hubble 
for a. Thus, an a priori fixing of w c would result to an 
inconsistent system [49]. 

To study the cosmological evolution, wc look at the 
system 



-3(1 + w c )p m H 



(5.24) 
(5.25) 



where H — a /a is the Hubble parameter. We will analyze 
the above system of the isotropic case for Kb — and 
Ab ^ 0, because of the different features that arise in the 
two choices of this parameter. 



5-2-1. Evolution of the System for = 

From the above dynamical system, taking into account 
(5.22), we find that there is only one fixed point in the 
evolution, the one with 



( Pm ,H 2 ,w c ) = (0,0,1/3) . 



(5.26) 



Linear perturbation around this point reveals that it is 
an attractor. From (5.23) and the conservation equation 
(5.25) wc find for the Hubble parameter 



a 
a 



2w r 



(l + 3w c )(l-3w c )(l + w c ) ' 



(5.27) 



The above equation can easily be integrated and solved 
for w c (the "— " sign in the solution of w c is rejected 
because it gives imaginary Hubble parameter) 



1 2 a 4 
Wc_ ~3 + 3 VfT^ 



(5.28) 



From this equation, we see that from any initial con- 
dition along the line of tuning (5.22), the expansion of 
the universe drives the equation of state to w c — > 1/3, 
i. eradiation. We have also verified this by integrating 
the system numerically. During this cosmological evolu- 
tion, it can never happen that w c = wf (compare (5.21) 
with (5.22)) and thus the whole system is regular. 



5-2-2. Evolution of the System for Kb / 

From the dynamical system, as written in the previ- 
ous subsection, taking into account (5.22) with Kb ^ 0, 
we find that there is a fixed point in the cosmological 
evolution 



{p mi H 2 lWc ) = p f 



Pf 



-1 



(5.29) 



with 



a Pf _ 
m£7 



1+3 %ff- ■ Since we should have 
this fixed point exists only for Kb > and 



Pm > 

corresponds to a de Sitter vacuum. Linear perturbation 
around this point reveals that it is an attractor. Thus, 
any matter density on the brane eventually evolves to a 
state of vacuum energy. 

A potentially interesting case, with a cosmological evo- 
lution resembling that of our universe, would be the one 
where < aKs/M^ <C 1. Then, the line of isotropic tun- 
ing, as it is given by relation (5.22), has a maximum 
close to w c ~ 1/3. The evolution of an initial energy 
density larger than the one corresponding to that max- 
imum, will evolve towards w c <~ 1/3, pass from w c = 
and asymptote to w c ^ -1 (see Fig. 1 for an example). 
The asymptotic value of the effective cosmological con- 
stant at the fixed point will be a^M- ~ <C 1 and 

PI M 6 

should be rather small to match with observations (when 
performing such a comparison, it is reasonable to assume 
that all the dimensionful scales of the theory are roughly 
of the same order, i.e.a^ 1 ^ 2 <~ Mpi ~ Mq). The lat- 
ter requirement of extremely small aAg /Mq is the usual 
cosmological constant problem. Although the standard 
cosmological evolution is described by piecewise constant 
equations of state with w ~ 1/3, 0,-1, in the present the- 
ory the equation of state has always time dependence. 

For Kb < 0, there does not exist any fixed point. The 
evolution of this system has a runaway behaviour and 
flows to w c — > oo while p m — > + . During this evolution, 



the singular point 



can be encountered only 



if K B I < -.326 (equating the expressions (5.21) and 
(5.22) the energy density p m can be real and positive only 
for this range of Kb). It is easily verified that for these 
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0.5 



-0.5 




400 



FIG. 1: The evolution of the equation of state w c for 
the case of aA^/Mg = 1CT 4 with initial conditions 



31 and ffi = .02. 



pi 



"dangerous" values of As, it is > 1/3. Therefore even 
in the case in which Ab takes these values, the dangerous 
point is reached only if initially w c < io+, 

We finally note that, for As ^ 0, there does not exist 
a fixed point with (p m ,H 2 ) = (0, 0) because the equation 
of state w c diverges at that point. 



6. THE UNCONSTRAINED 
ANISOTROPIC CASE 

In the previous section we found some interesting cos- 
mological evolutions of an isotropic brane-universe if 
Ab > 0. Wc should keep in mind however, that in all 
cases studied, the energy density on the brane is tuned to 
the equation of state in a specific way, which seems at first 
sight artificial. Therefore, it would be worth studying 
some cosmological evolution, in a codimension-2 branc- 
world model, in which this tuning is not required. If the 
system then evolves towards the previously studied line 
of isotropic tuning, we will conclude that this tuning is 
an attractor and thus not artificial. 

To do so, we have to consider geometries which are 
not isotropic and in which the Riemann tensor cannot 
be expressed in terms of the Ricci tensor and the curva- 
ture scalar [48]. Then, the constraint equation (5.8) will 
not give rise to a brane-bulk matter tuning, but rather 
to a dynamical equation for the anisotropy of the space. 
For this purpose, let us consider the following anisotropic 
ansatz where the metric functions depend only on the 
time t and the radial coordinate r (i.e. keeping the az- 
imuthal symmetry) 

3 

ds 2 = -N 2 (t,r)dt 2 +^2A 2 l {t 7 r){dx i f+dr 2 +L 2 (x,r)d6 2 . 
i=l 

(6.1) 

We can again use the gauge freedom to fix N(t, 0) = 1, 



while we define Ai(t,0) = a,i(t). The singularity con- 
ditions dictate as before that N'(t,0) = A<(*,0) = 0, 
while the second derivatives of these metric functions are 
unconstrained. The most general anisotropic evolution 
scale factors with the above property can be written as 

A 1 (t, r) = a(t)b(t)c(t) + (t)r 2 + ... , (6.2) 



A 2 (t,r) = ^+^(ty + ... 

A 3 (t,r) = ^+Ut)r 2 + ... 
c{t) 



(6.3) 
(6.4) 



where a = (aia 2 a 3 ) 1 / 3 represents the "mean" scale fac- 
tor and b, c represent two degrees of anisotropy. To sim- 
plify further the analysis of the system, we will choose 
c = const, (by a coordinate redefinition then we can al- 
ways set c = 1). The dynamics of this particular choice 
can help us to understand the qualitative features of the 
general case. 

The purely four-dimensional brane Einstein equations 
take the form 



a 2 b : 
3-2 -- 
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(6.5) 
(6.6) 
(6.7) 



while the equation which is coming from the extra di- 
mensions is now dynamical, providing a Hubble equation 
for b 



b 2 
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(6.8) 



We will now assume, as before, that the vacuum en- 
ergy (tension) of the brane cancels the contribution A4 
induced my the deficit angle as in(5.16). After this sim- 
plification, the above equations read 



a 2 b 2 M 2 pl ' 

o d , a" 2 ft \ P m 
a a 2 M 2 
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b_tf ab 
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(6.9) 
(6.10) 
(6.11) 



while the equation coming from the extra dimensions be- 
comes 



b 2 



4M 2 P[ 



3 
32q^ 
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= f(Pm,w) 



{3w - 1) + 2(2w + l)a 



(6.12) 



It is interesting to observe that the Hubble equation (6.9) 
for the "mean" scale factor a, after substitution of (6.12), 
has apart from the usual linear term in p (of the conven- 
tional four-dimensional cosmology) , additional correction 
terms in p. This is similar to what happens also to five- 
dimensional brane-world models [9] and is due to the 
presence of extra dimensions. This modification hap- 
pens only in the anisotropic case. In the pure isotropic 
case the four-dimensional brane-universe feels the extra 
dimensions by only adjusting its energy density to its 
equation of state, but without any modification in the 
structure of the Fricdmann equation. 



6-1. Cosmological Evolution of an 
Anisotropic Brane-Universe 

Before analyzing the system, let us note again that 
in contrast to pure four dimensional anisotropic dynam- 
ics, where a constant w is allowed, in the present case, 
where there is an extra dynamical equation because of 
the presence of the bulk, w has to evolve. Indeed, by 
differentiating (6.12) and then using the four dimen- 
sional equations of motion and the conservation equation 
p m + 3(1 + w)p m ^ = we can find a differential equation 
for w 



dw 



w + 3 



2/ - (1 + w)p n 



df 



dp„ 



= . 



(6.13) 



Imposing a further condition to keep w constant, would 
result to a constant p m related to w in a specific way and 
by the conservation equation, to zero Hubble for a. Thus, 
an a priori fixing of w would result to an inconsistent 
system [38, 49]. 

To have real Hubble parameters for a and b, p m and 
w have to lie in specific regions for which the following 
inequalities are satisfied 



/>0 



/ + 



Pr, 



M 2 Pl 



> . 



(6.14) 



Define the following boundaries of the allowed regions 



-2 



+ 



Wl = 



i 1 - 



U>2 = 



3^(1 + !^) 



(6.15) 



(6.16) 



The quantity w\ coincides with the line of isotropic tun- 
ing w c given in (5.22). The inequalities (6.14) can then 
be re-expressed as conditions for p m and w: 



• For p m > 0, we should have w > w±. 

• For p m < 0, we distinguish two cases 

- for — | < jjfir < 0, we should have w < W2- 

- for j£*r- < — |, we should have w > w 2 . 

To study the cosmological evolution, we look at the 
system 



2 V ' M 2 



f , 



pi 



Hb — —3H a Hb , 

Pm = -3(1 + w)p m H a , 



(6.17) 

(6.18) 
(6.19) 



where H a = a/a and Hb — b/b are the Hubble parameters 
for a and b respectively. The third equation is the energy 
conservation equation in which only the Hubble param- 
eter H a for the "mean" scale factor a appears. Again we 
will analyze the anisotropic case for = and Ag ^ 
and we will compare the cosmological evolution with the 
cosmological evolution of the tuned isotropic case. 



6-1-1. Evolution of the System with Kb = 

There are three different regions in the (w,p m ) plane 
where these inequalities are satisfied. This relative free- 
dom to choose the matter on the brane should be com- 
pared to the tuning that happens in the isotropic case. 
Relaxing the isotropy condition, the system can have ini- 
tial conditions in a vast region of the parametric space. 

From the dynamical system (6.17), (6.18), (6.19) we 
see that there is only one fixed point and that it is the 
same with that of the isotropic evolution, i.e. 



(p m ,H 2 a ,Hlw) = (0,0,0,1/3) 



(6.20) 



Linear perturbation around this point reveals again that 
it is an attractor. 

The presence of the previous attractor fixed point 
will drive the system towards a final isotropic state of 
radiation. The way in which this fixed point is ap- 
proached from an arbitrary initial energy density, can 
tell us whether the line of isotropic tuning is an attractor 
or not. If the anisotropy monotonically decreases during 
the evolution, it means that the line of isotropic tuning 
is an attractor. 

In order to analyze the features of the anisotropic evo- 
lutions, we proceed numerically. We solve the system of 
the two second order equations (6.10), (6.11) and the two 
first order equations (6.13), (6.19) for the four functions 
a, b, p m , w. To understand how the anisotropy involves 
we define the mean anisotropy by the following quantity 



\l 
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((H) -Htf 
3(H) 2 



(6.21) 
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where Hi — dj /en (with m defined after (6.1)) and (H) = 

Our analysis shows that, the relaxation of the tuning 
relation between w and p m , has as a consequence the ap- 
pearance of new branches of brane world evolution, while 
the system tends quickly to the isotropic fixed point at- 
tractor. Furthermore, we observe that the anisotropy in 
all cases decreases much more quickly than in the four di- 
mensional case with the same initial conditions but with- 
out the extra dimensional constraint. 



6-1-2. Evolution of the System with Ab / 

Studying the asymptotics for w\ and u> 2 from (6.15), 
(6.16) we find that for Ab ^ 0, there are three intervals 
of Ab with different shape of the allowed regions in the 
{p m ,w) plane. 

From the dynamical system (6.17), (6.18), (6.19) we 
see that there are two fixed points in general. The first 
one is the same with that of the isotropic evolution, i.e. 

(p m ,H 2 a ,H 2 ,w) = (^,^,0,-1) , (6.22) 

with j£f = -f + fy/l + I ^ and it exists only for 

Ab > 0. Linear perturbation around this point reveals 
again that it is an attractor. 

The second fixed point that we find is 

( Pm ,HlHlw)= (-A S ||, 0,||,l) • (6-23) 

and it exists only for A B > 0. Linear perturbation 
around this point reveals that it is a repeler. 

Whenever the previous attractor fixed point exists, 
the system will be driven towards a final isotropic dc 
Sitter state. The way in which this fixed point is ap- 
proached from an arbitrary initial energy density, can 
tell us whether the line of isotropic tuning is an attractor 
or not. We analyzed the anisotropic system numerically 
and we found [49] that, the relaxation of the tuning be- 
tween w and p m , has as a consequence the appearance of 
new branches of brane world evolution, while for choice 
of the parameters the system tends to the attractor fixed 
points, whenever they exist. For the cases when Ab > 
the lines of isotropic tuning are attractors, with strength 
depending on the value of A^. 

Furthermore, we observe that the anisotropy in all 
cases, apart for Ab < 0, decreases much more quickly 
than in the four dimensional case, with the same initial 
conditions, but without the extra dimensional constraint. 
For Ab < 0, the anisotropy increases and is larger than 
the one of the purely four dimensional case. 

Let us now examine again the interesting possibility of 
< ccAb/Mq <C 1 with the transition between w = 1/3 
to w = and finally to to = -1. As it can be inferred 
from the previous discussion, the system tracks the line of 



isotropic tuning and evolves towards the attractor fixed 
point with w = — 1. However, due to the small value 
of aAs/Mg, the line of isotropic tuning is a very weak 
attractor. Most of the evolution is rather anisotropic with 
A <~ O(l) — O(10 _1 ) until the fixed point is approached, 
in which region it drops to zero. 

This large anisotropy makes the cosmological evolu- 
tion phenomenologically problematic. In order that the 
anisotropy is acceptably small, the initial conditions for 
the energy density and the equation of state should be 
fine tuned to lie very close to the line of isotropic tuning 
initially. 

In conclusion, by analyzing the anisotropic dynamics of 
the system we have seen that the lines of isotropic tuning 
are attractors for Ab > 0, with A^-dependent attracting 
strength. The most phenomenologically accepted evolu- 
tions with < oAs/Mg < 1 do not isotropise quickly 
enough and thus need a fine tuning in order to evolve 
with acceptable anisotropy. 



7. SIX-DIMENSIONAL 
SUPERGRAVITY COSMOLOGICAL 
MODELS 

There are also cosmological models in six dimensions 
coming from supergravity theories [53]. In some early 
works cosmological solutions were found in D=6, N=2 
supergravity theories. These solutions describe vac- 
uum state or radiation four-dimensional universe in a 
Minkowski xS 2 space [42]. 

Chaotic inflation was studied in an anolaly-free, 
gauged (1,0) supergravity model in six dimensions in 
which the unique vacuum state is Minkowski space- 
time cross an internal S 2 which leaves half of the six- 
dimensional supersymmetries unbroken [54]. It was 
shown that inflationary dynamics consistent with the cos- 
mological constraints can be realized provided that the 
radius of the internal space S 2 satisfies a constraint. In 
this model, the inflaton field <j> originates from the com- 
plex scalar fields in the D=6 scalar hyper- multiplet. The 
mass and the self couplings of the scalar field are dictated 
by the D=6 Lagrangian. The scalar potential has an ab- 
solute minimum at <\> = with no undetermined moduli 
fields. 

The model is based in a Lagrangian which have been 
constructed in [55, 56]. It is chiral, hence potentially 
inconsistent due to the presence of gauge, gravity and 
mixed anomalies. In [41], an anomaly free model with 
a gauge group E§ x E-j x U{1)r has been constructed 
which involves a hyper-scalar multiplet transforming in 
a 912-dimensional pseudo-real representation of £7. In a 
convenient parameterization, the potential for the scalars 
take a simple form 

V = 4^H(1 + \4>?f + 9 \{4>T a m (7.1) 
K Si 
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where n is the D=6 gravitational coupling and has a di- 
mension of square length. The g\ and are the cou- 
pling constants of the U(1)r and £7, respectively, and 
they have the dimensions of a length, the T a are the 
Hcrmitian generators of £7 in its 912-dimensional repre- 
sentation. The most important property of this potential 
is that it has a unique minimum at (f> = 0. Chaotic in- 
flation is realized for the zero KK mode with a potential 
(7.1) without the constant term. 



8. CONCLUSIONS 

In this talk we discussed cosmological models in six 
dimensions. First we presented a (4+l)-braneworld cos- 
mological model in a six-dimensional bulk. If a(t) = b(t), 
with a(t) the usual scale factor of the three physical di- 
mensions, and b(t) the scale factor of the extra fourth di- 
mension, we found the generalized Friedmann equation in 
six-dimensions of the Randall-Sundrum model describing 
the cosmological evolution of a four-dimensional brane- 
universe. If a(t) 7^ b(t) the four-dimensional universe 
evolves with two scale factors. However, for an observer 
in the moving brane, a and b are static depending only 
on the coordinate on which the 4-brane is moving. Then, 
demanding to have an effective Friedmann-likc equation 
on the brane, we showed that the motion of the 4-brane 
in the static bulk is constrained by Darmois-Israel bound- 
ary conditions, resulting in a relation connecting a and b 
acting as a constraint of the brane motion. 

The way a and b are related depends on the energy- 
matter content of the 4-brane. We then explored what 
are the consequences of the presence of this constraint 
of motion for the cosmological evolution of the branc- 
universe. We assumed that the universe started higher- 
dimensional at the Planck scale with all the dimensions 
at the Planck length, and subsequently an anisotropy was 
developed between the three physical dimensions and the 
extra-dimension. We then followed the evolution by mak- 
ing a "phenomenological" analysis of how the two scale 
factors evolved under various physical assumptions. In all 
considered cases, (A)dS and Minkowski six-dimensional 
bulk, open, closed and flat brane-universe and matter, 
radiation, cosmological constant and dark energy domi- 
nated three-dimensional physical universe, we found that 
dark energy is needed for the dynamical suppression and 
subsequent freezing out of the extra fourth dimension. 

We next discussed the gravitational dynamics of con- 
ical codimension-2 branes of infinitesimal thickness. We 
have considered theories with codimension-2 branes 
which arc augmented with an induced gravity term on the 
brane and a Gauss-Bonnet term in the bulk. We made 
the observation that the dynamics of the these theories 
is not exhausted by studying the "boundary" Einstein 
equation, which is exactly four-dimensional and bears no 
information of the internal space (modulo a cosmologi- 
cal constant contribution). In the case of a pure brane 
induced gravity term, the higher dimensional Einstein 



equations evaluated at the position of the brane, give a 
very precise and strong relation between the matter on 
the brane and the matter in the bulk in the vicinity of 
the brane. In other words, the bulk energy content is 
the primary factor for the cosmological evolution on the 
brane. Alternatively, for a static matter distribution on 
the brane to be possible, there should exist its bulk mat- 
ter "image". 

This strong relation, that we have noted, can be 
avoided with the inclusion of a bulk Gauss-Bonnet term. 
An even more natural way to achieve this is to relax 
the requirement of purely conical branes and admit gen- 
eral brane solutions with an appropriate regularization 
(thickening of the brane), so that the singularities are 
smoothened. In view of the difficulties related to the 
Gauss-Bonnet term in the context of effective field the- 
ory, we pointed out that the thickening of the brane is the 
most physical direction that one should follow in order 
to discuss the dynamics of codimension-2 branes. 

We then discussed the cosmological dynamics of a con- 
ical codimension-2 braneworlds in a theory with a Gauss- 
Bonnet term in the six-dimensional bulk and an induced 
gravity term on the three-brane. For simplicity, we con- 
sidered that the bulk matter consists only of a cosmo- 
logical constant As but the brane matter is general and 
isotropic. We then analyzed in detail the Einstein equa- 
tions evaluated on the boundary. 

We studied the system first for an isotropic metric 
ansatz. In the pure induced gravity dynamics there is 
a tuning between the matter allowed on the brane and 
in the bulk. This tuning, when the matter in the bulk 
is only a cosmological constant, gives a precise relation 
between the matter density on the brane and its equa- 
tion of state. If a Gauss-Bonnet term is added in the 
bulk, the constraint equation giving the previous tuning 
is modified by the addition of a Ricmann squared term. 
However, since for isotropic evolutions the Ricmann ten- 
sor can be expressed in terms of the Ricci tensor and 
the scalar curvature, the conclusion about the presence 
of the tuning remains the same as in the induced gravity 
case. We found that if As > the system has a fixed 
point with equation of state w = — 1 and corresponds to 
a de Sitter vacuum. If As =0, the system has a fixed 
point with equation of state w = 1/3. Both of the pre- 
vious fixed points are attractors. If As < the system 
has no fixed point and the evolutions exhibit a runaway 
behaviour to w — > +00. 

We then looked on the cosmological evolution of the 
system for a particular anisotropy. If the system starts its 
evolution in the region of parametric space which has as 
a boundary the line of isotropic tuning, it tracks the lat- 
ter line and isotropises towards the attractor fixed point 
with w = —1 for Kb > 0, or the one with w = 1/3 for 
Ab = 0. In the two other regions the system has a run- 
away behaviour w — > 00 and does not isotropisc. For 
Ab < 0, the system has always a runaway behaviour. 
The important result of this analysis is that the line of 
isotropic tuning, is an attractor. However, for values of 
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Ab which give acceptable cosmological evolutions, a fine 
tuning is unavoidable because of the weak strength of the 
above-mentioned attractor. 

Six-dimensional supergravity models have a much 
richer structure. There are many scalar fields with a 
scalar potential fixed by the supergravity theory. There 
are also gauge fluxes which guarantee the stability of 
the theory. Recently, they are six-dimensional warped 
braneworld solutions in supergravity teories [57] and 
it would be interesing to find also codimension-1 or 
codimcnsion-2 cosmological braneworld solutions in six- 



dimensional supergravity theories. 
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